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Abstract 

We study the bond percolation problem under the constraint that the total number of occupied 
bonds is fixed, so that the canonical ensemble applies. We show via an analytical approach that at 
criticality, the constraint can induce new finite-size corrections with exponent j/can = '^yt — d both 
in energy-like and magnetic quantities, where yt = l/i' is the thermal renormalization exponent 
and d is the spatial dimension. Furthermore, we find that while most of universal parameters 
remain unchanged, some universal amplitudes, like the excess cluster number, can be modified and 
become non-universal. We confirm these predictions by extensive Monte Carlo simulations of the 
two-dimensional percolation problem which has ?/can = 
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I. INTRODUCTION 



model 



The bond-percolation model [l[ can be considered as the g — )■ 1 limit of the g-state Potts 



3[. Consider a lattice G = iV^E) with (.E) the vertex (edge) set, the reduced 
(divided by kT) Hamiltonian of the Potts model reads 

'H{K,q) = -KY,^<^.^, (a = l,---,g), (1) 

where K is the coupling strength. By introducing bond variables on the edges and integrating 
out the spin degrees of freedom, one can map the Potts model onto the random-cluster (RC) 



model. This is the celebrated Kasteleyn-Fortuin transformation j4|. The partition sum of 
the RC model thus assumes the form 

2rc(n,g) = E ^'^'^'^-"^ (n = e^-l), (2) 

where the sum is over all the spanning subgraphs A of G, |^| is the number of occupied 
bonds in A, and k{A) is the number of connected components (clusters). The Kasteleyn- 



Fortuin mapping provides a generalization of the Potts model to non-integer va! 



ues g > 0. 



In the g — )■ 1 limit, the RC model reduces to the bond percolation problem . In this 
limit, the partition sum assumes the non-singular form {u + 1)'^', where \E\ is the number 
of the lattice edges. Nevertheless, rich physics still exists, e.g. in the derivatives of the 
partition sum involving g. 

Percolation in two dimensions has been extensively studied. Percolation thresholds on 
many two-dimensional lattices are exactly known or have been determined to very high 
precision (see e.g. Ref. jsj] and references therein). Most of the critical exponents are also 
exactly known. For instance, the leading and subleading thermal renormalization exponents 
are yt = l/v = 3/4 and yi = —2, respectively; the leading two critical exponents in the 



There still exist some 



magnetic sector are yh = d — j3 / u = 91/48 and yh2 = 71/96 
critical exponents whose exact values are unknown, including the backbone exponent and 
the shortest-path fractal dimension Much progress has recently been made in the 



context of the stochastic-Lowner evolution 
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111. 



The present work provides a study of another aspect of percolation. We introduce an 
"energy-like" constraint that the total number of occupied bonds is fixed and study the 
effects of such a constraint on the critical finite-size scaling (FSS). Since occupied bonds can 



formally be considered as particles, we shall refer to percolation under this constraint as 
percolation in the canonical ensemble, and to the case without the constraint as percolation 
in the grand-canonical ensemble. Two of us have, for several years, been studying the 



effects of such constraints 
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-|l4j. The leading FSS in systems under energy-like constraints 



is derived on the basis of the so-called Fisher renormalization [15| procedure, which was 
originally formulated for systems in the thermodynamic limit. Consider an energy- density- 
like observable e in the critical RC model that scales as {e{L)) = + EsL^*''^ in the grand- 
canonical ensemble, where L is the linear system size and d is the spatial dimensionality. 
The constraint in the canonical ensemble implies: 

(i) , e{L) for a given size L is fixed at the expectation value Ea in the thermodynamic limit. 
This expectation value is normally different from the grand-canonical finite-size average 
{e{L)). 

(ii) , fluctuations of e are forbidden — i.e., (e^) = {e)^. 

The effect of (i) is accounted for by the Fisher renormalization. The bond density l^l/li?! 
is such an observable. For bond percolation, however, since the bond variables on different 
edges are independent of each other, = {\A\)/\E\ does not depend on the system size 
L — i.e., if we write pb = Pb,a + Pb,sL^*~'^, then the amplitude pb^s = 0. Thus, percolation 
provides an ideal system to study the effects of energy-like constraints due to the suppression 
of fluctuations. 

The remainder of this work is organized as follows. Section II describes the sampled 
quantities and the FSS of physical observables for percolation in the grand-canonical ensem- 
ble. Section III derives the effects of the "energy-like" constraint in the critical FSS. The 
numerical results are presented in Sec. IIVI A brief discussion is given in Sec. |Vl 



II. SAMPLED QUANTITIES AND FINITE-SIZE SCALING 
A. Simulation and sampled quantities 

We consider the bond-percolation problem on an L x L square lattice with periodic 
boundary conditions. The grand-canonical simulations follow the standard procedure: each 
edge is occupied by a bond with probability p, after which the percolation clusters are 



constructed. For simulations in the canonical ensemble, a Kawasaki-like scheme 16| is used. 
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Given an initial configuration with a total number of occupied bonds |^| = p\E\, where 
lE'l = 2L^, an update is defined as the random selection of two edges and the exchange of 
their occupation states. Each sweep consists of 2L^ updates. Quantities are sampled after 
every sweep. Note that although the bond number |^| = p\E\ happens to be an integer for 
the critical square-lattice bond percolation, p\E\ is generally a real number. In this case, 
one can simulate at the ceiling and the floor integer of pl-El and apply an interpolation. 

Given a configuration A as it occurs in the simulation, we denote the sizes of clusters as 
Ci {i = 1, . . . , k{A)); Ci is reserved for the largest cluster. The following observables were 
sampled. 

1. Energy-like quantities: 

• The bond-occupation density pf, = 

• The cluster-number density pk = {k{A)) /\V\. 

2. Specific-heat-like quantities: 

. a = L~'^{{\A\')-{\A\y). 
. C, = L~''m^?)-{kiA))'). 

. C2 = L-^ {{[\A\/2 + kiA)f) - {\A\/2 + k{A)f) =Ck + C^t + a/4. 

Here, the factor 1/2 in the definition of C2 arises from the fact that the critical line 
Uc{q) of the RC model ([2]) on the square lattice is described by Uc{q) = g^^^, which 
leads to a weight gl-^l/^+^^C-^) for a subgraph A along the critical line. In the canonical 
ensemble, the bond number |^| is fixed, and Cb and Ckb reduce to 0. 

3. Magnetic quantities: 

• The largest-cluster size Si = (Ci). 

• The cluster-size moments Si {i > 1). Defining Sg = Cf , we sampled 5*2 = (1S2) 
and (35I-254). 

4. Dimensionless quantities: 
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• Wrapping probabilities Rb, Rx, and Re. The probability Ri counts the events 
that configuration A connects to itself along the x direction, but not along the 
y direction; Rf, is for simultaneous wrapping in both directions; R^ is for the 
X direction, irrespective of the y direction; R^ is for wrapping in at least one 
direction. They are related as Re = 2Ri + Ri, and R^ = Ri + Rb, thus only two 
of them are independent jl^, 

• Universal ratio Qs = {Cf)/Si. 

• Universal ratio Qm = {^2)'^/ (35| — 2S4). 



B. Finite-size scaling in the grand-canonical ensemble 

The critical FSS of the sampled quantities can be obtained from derivatives of the free- 
energy density / = —L~'^\nZy.c (^rc is given by Eq. ([2])) with respect to the thermal scaling 
field t, the magnetic scaling field h, or the parameter q. In the grand-canonical ensemble, 
the FSS of /(g, t, h, L) is expected to behave as 

/(g, t, h, L) = fr{q, t, h) + L-^fs{q. tLy% hLy\ 1) , (3) 

where higher-order scaling fields have been neglected, and fr and fs denote the regular 
and the singular part of the free-energy density, respectively. The thermal scaling field t is 
approximately proportional to u — Uc in Eq. ([2]), where Ue{q) = y/q is the critical line of the 
g-state RC model on the square lattice. 

FSS for energy-like quantities. From the partition sum in Eq. ([2]) and Eq. Q, it can be 
derived that at criticality 

(-g)(d/,/dg) = L-''{\A\/2 + k{A)) = pb{L) + pu{L) = pb,o + Pk,o + bL-"" , (4) 
{-u/2){df/du) = 2-'L-''{\A\) = pb{L) = pb,o + aiy^-'' , (5) 
{-q){df/dq) = L~\k{A)) = pkiL) = pk,o - aL^^"^ + feL"'^ , (6) 

where fc = f{u = Mc) is the free-energy density along the critical line Uc{q) = y/q. The 
last equality in Eq. (jlj) refiects the analyticity of fc (including the amplitude of its finite-size 



dependence) along the critical line 19| . Moreover, the correction amplitude b is universal 19 



20], which is also referred to as the excess cluster number for percolation. The term p^ Q 



5 



follows from the exact results for the critical Potts free energy on the square lattice 21|. 



22|. The last equality in Eq. 



An exact result is also available for the triangular lattice 
arises from Eq. ([3]), and pf,o accounts for the background contribution; the self-duality of 
the square-lattice model yields Pb,o{(l) = 1/2- In the g — ?■ 1 limit, the amplitude a vanishes 
as a ~ g — 1. Therefore, one has for critical percolation 

Pb{L) = pbfi , pk{L) = pkfl + bL''^ . (7) 

FSS for specific-heat-like quantities. Similarly, one has the second derivatives (at q = I) 



as 



d^^ 



q^] fc = C2 = C2,o + cL-'^ , (8) 



u-^] f = a = a,o , (9) 



d \ f d 

du 



(^q^^ f = Ckb = Cktfi + 2a' ly^-^ + c'L-'' , (10) 
- (^q^^ f = Ck = Ck,o - 2a'L^'-'^ + c'L"'^ , (11) 
where the last equalities in Eq. f lTOj) and f|TT]) arise from the fact that the parameter a in 



19 



is 



19|, 



Eqs. (|5]) and ([6]) behaves as a ~ a'{q — 1). The analyticity of fc along the critical line 
also reflected in Eq. (|8]). It is also derived that the amplitude c is a universal quantity 
called the excess fluctuation of the number \A\/2 + k{A) for critical bond percolation. 

FSS for magnetic quantities. The critical FSS of magnetic quantities can be obtained by 
differentiating the free-energy density with respect to the magnetic scaling field h. It can be 
shown that at criticality, 

Si{L) ~ L^^ , S2{L) ~ L^y^ . (12) 

Corrections to FSS. At criticality, the asymptotic behavior of a quantity 0{L) is supposed 
to follow the form 

0{L) = L^{0o + corrections), (13) 

where ip is the leading critical exponent and Oq is the amplitude. For our observables in 
bond percolation, values of ip are given by the previous analysis {tp = for dimensionless 
quantities). The FSS theory predicts several types of "corrections" in Eq. (fT3!) (see e.g. 



Refs. 



231] and 2J] for a review and references therein). These terms include corrections 
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from the irrelevant scahng fields; the leading one has exponent Ui = —2 for percolation. A 
regular background term, which appears e.g. as L'^^ for 5*1 and L*^"^^*" for 5*2, contributes 
to the "corrections". These regular background terms are also present in the associated 
dimensionless ratios. For the energy-like quantities (p^ and pb) and the specific-heat-like 
quantities, ip is negative, thus the regular background (which is L'"^ here) included in 
the "corrections" in Eq. (fT3!) is actually not a correction, instead it describes the leading 
behavior, which is the thermodynamic limit of the sampled quantity. 



III. FINITE-SIZE SCALING IN THE CANONICAL ENSEMBLE 

As mentioned in the Introduction, the critical FSS for bond percolation in the canonical 
ensemble is due to the suppression of the fiuctuation of occupied-bond density pb- In this 
section, we demonstrate how the critical FSS is affected by such a suppression. 

Without the constraint, any configuration A with bond number |^| = A*";, occurs with 
probability p^''(l — p)(^<=~^b) (^N^ = \E\). Thus, the value Og of an observable O in the 
grand-canonical ensemble can be written as 

0,(p,L) = 5^p^^(l-p)^-^^ ^ 0{A) 

Nt=0 A:\A\=Nb 

(14, 

where p = Nb/N^ is the occupied-bond density and Oc is obtained by averaging over all the 
configurations A with \A\ = Nb, of which the total number is NJ/Nb\{Nf. — Nb)\. Namely, Oc 
is the canonical-ensemble value of observable O. For a sufficiently large system, the binomial 
distribution in Eq. (fT4l) is well approximated by a Gaussian distribution 

/(p,p,L) = ^e-(^-^)V2-\ a=^pil-p)/Nc, (15) 
Vivrcr 

where p = p is the average bond density for occupation probability p, and the variance a 
decreases as oc 1/\^N^ oc L"'^/^. Thus, Eq. f lT^ can be re- written as 



Og{p,L) = / dpf{p,p,L)Oc{p,L) 



(16) 



In principle, based on the known scaling behavior of Og in the grand-canonical ensemble, 
the critical FSS of Oc can be obtained by the inverse orthogonal transformation of Eq. f|T6|) . 
Here, we take an approximation approach by Taylor expanding Oc(p, L) in Eq. (ITB]) in powers 
of 6p = p — J) about p = p = p, and term-by-term evaluation of the integrals: 

J dp f{p, p, L)OM L) = 0,{p, L) + 0',{5p)f + 0:{{5pf/2)f + ... (17) 

where the derivatives O'^ and O'^ are taken at p, and the average ( )/ is over the Gaussian 
distribution in Eq. ( TT5l) . It is easily derived that {5p) j = and ((5p)^/2)/ = p(l — p)/2Ne. 
Combination of Eqs. ( IT6l) and (ITTl) yields 

Og{p,L)c^O,{p,L) + BO';L-'' , (18) 

with i? = p{l — p)/4 a non-universal constant. 

Near criticality pc, we denote Ap = p — Pc- Since p = p, we can equivalently denote the 
distance to the percolation threshold as Ap = p — Pc with p^ = pc and Ap = Ap. Thus 

Og{Ap, L) = Og,r{Ap) + L^[Og,s{Ly^Ap) + 0,,L^»], (19) 
Oe(Ap, L) = 0,,,(Ap) + L^' [©...(L'^" Ap) + ^eani^^^^" + O^^iiy^] , (20) 

where new finite-size corrections with exponent i/can are allowed in Eq. (!20l) . It is known 
from Eq. ([20]) that 0;'(p„ L) = 0^^ + L'^'+'^y^Ol^. Taylor-expanding the r.h.s of Eqs. dUEO]) 
and substituting them in Eq. (fT8|) give 

O,,,(0) + L^Pg^M + O'^^^iy^Ap + O.^L^^' + ...]= (21) 

o,,,(o) + L'^'[o,,,(o) + o;,L^''Ap + Oe,ii:^' + BO'i^i^y"-'' + 5,ani^^^^" + bo'I.l-''-^' + ...] 

Assuming that 2yp — d < 0, the term with O"^ acts as a finite-size correction. By 
comparing the scaling behavior of the l.h.s and the r.h.s. of Eq. (^T^, we obtain 

1. tjj' = ^1; and Vp = yt, 

2. Oc,r{0) = Og^riO), Oc,s(0) = Og,,(0) (for ip 7^ -d), and ^ = O'^ ^ etc. In other words, 
the scaling functions in Eqs. ( 1T9|) and ( l20l) are identical up to some correction terms. 

3. new correction terms appear in the canonical ensemble, and they have exponents 
2yt — d and —d — ip (if —d — ip < 0). This is because the leading finite-size correction 
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in the l.h.s of Eq. ( 12T|) is described by an exponent Ui = —2, and thus the terms with 
l^vp-d _ ijiyt-d j^-d-^p' _ j^-d-i(! i^oyQ cancelled by the newly included 

correction terms with It can be shown that corrections with exponent n{2yt — d) 

with n = 2,3,... can also occur if higher-order terms are kept in the above Taylor 
expansions. 

4. for the case of ip = —d, Oc,s(0) = Og,s(0) does not hold. Instead, one has Og,s(0) = 
Oc,s(0) + BO" J.. This predicts that the excess cluster number b in Eq. ([7]) is changed 
by the constraint; further, since B and p"^ are non-universal, the canonical-ensemble 
value of b is no longer universal. 

Note that the assumption 2yp — d < holds since 2yt — d < for percolation in two and 
higher dimensions. 

Finally, we mention that although the constraint is "energy-like" , the derivation applies 
to both the energy-like and the magnetic observables. 

IV. NUMERICAL RESULTS 

To examine the above theoretical derivations, we carry out Monte Carlo simulations for 
the critical bond percolation on the square lattice, with a fixed bond number = Pc\E\ = 
L^, and for 21 sizes in range 4 < L < 4000. The number of samples was about 10^ for 
L < 480, 5 X 10^ for L = 800, 2.5 x 10^ for L = 1600, and 10^ for L = 4000. For the purpose 
of comparison, additional simulations were also carried out in the grand-canonical ensemble. 

A. Evidence for the correction exponent yean = 

2yt-d 

We first examine the critical FSS of the wrapping probabilities. In the grand-canonical 
ensemble, the finite-size corrections arise only from the irrelevant scaling fields for which 
the leading exponent yi = —2. This is illustrated in the inset of Fig. [H In the canonical 
ensemble, the existence of the newly induced correction exponent ?/can = '^Ut — d = —1/2 is 
clearly demonstrated by the approximately linear behavior for large L in Fig. [H where -Ri 
is plotted versus 

According to the least-squares criterion, we fitted the data for these wrapping probabilities 
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FIG. 1: (Color online). Plot of the wrapping probability Ri versus L~^/^ in the canonical ensemble 
(main plot). The inset shows Ri versus in the grand-canonical ensemble. The straight lines 
represent the leading finite-size dependence as obtained from the fit. The error bars of the data 
are smaller than the size of the data points. 

by 

= 0o + Biiy^^-^ + SaL^^'^^" + Biiy' , (22) 

where Oq represents the universal value for L — oo, and the correction exponent Ui is 
fixed at —2. As a precaution against higher-order correction terms not included in the fit 
formula, the data points for small L < L^i^ were gradually excluded to see how the residual 

changes with respect to L^m- In general, we use results of the fit corresponding to a Lmm 
for which the quality of fit is reasonable, and for which subsequent increases of Lmin do not 
cause the value to drop by vastly more than one unit per degree of freedom. In practice, 
the word "reasonable" means here that is less than or close to the number of degrees of 
freedom, and that the fitted parameters become stable. 

The results are shown in Table [B The error margins are quoted as two times of the 
statistical errors in the fits, which also applies to other tables, in order to account for 



possible systematic errors. The universal values Oq can be exactly obtained 25|, l26|]. In 
the grand-canonical ensemble, these exact values of Oq were used and the amplitudes Bi 
and B2 were set at 0, so that the fit formula in Eq. fl22l) has only a single free parameter. 
Indeed, such a simple formula can well describe the data with L > 16 for all the wrapping 
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TABLE I: Fit results for wrapping probabilities. The exact values of Oq j25l. l26l|. and yi = —2 are 
used here. The abbreviations "GE" and "CE" stand for grand- and canonical ensemble, respec- 
tively. Entries "-" indicate the absence of a result, and the numbers without error bars are fixed 
in the fits. Error margins are quoted as two times of the statistical errors in the fits. 





Fit 


using = Oq + BiLy^'"^ 




+ Biiy^ 








Oo 


Bi 


yean 




B^ 


Vi 






GE 


0.169 415 435 


— 


— 


- 


0.23(1) 


-2 


16 






0.169 415 435 


— 


— 


- 


0.24(3) 


-2.02(4) 


8 


^1 




0.169 5(4) 


0.11(2) 


-0.53(5) 


- 


- 


- 


160 




CE 


0.169 415 435 


0.091 (8) 


-0.502(13) 


0.08(2) 


0.56(7) 


-2 


12 






0.169 415 435 


0.089 2(5) 


-1/2 


0.082 (4) 


0.55(3) 


-2 


12 




GE 


0.351642 855 


— 


- 


- 


-0.254(11) 


-2 


12 






0.351642 855 


— 


— 


- 


-0.30(5) 


-2.06(8) 


8 


Rb 




0.3515(6) 


-0.10(3) 


-0.55(7) 


- 


- 


- 


120 




CE 


0.351642 855 


-0.075(11) 


-0.50(2) 


-0.07(2) 


-0.53(5) 


-2 


8 






0.351642 855 


-0.075 8(8) 


-1/2 


-0.066(5) 


-0.53(3) 


-2 


8 




GE 


0.521058 290 


— 


— 


- 


-0.021 (16) 


-2 


16 






0.521058 290 


— 


- 


- 


-0.1(1) 


-2.4(6) 


8 


Rx 




0.5211(2) 


0.020 (4) 


-0.56(8) 


- 


- 


- 


32 






0.521058 290 


0.018(7) 


-0.54(8) 


- 


- 


- 


120 




CE 


0.521058 290 


0.014 (9) 


-0.51(9) 


0.016(13) 


- 


- 


16 






0.521058 290 


0.016(12) 


-0.52(11) 


0.01(3) 


0.03(8) 


-2 


10 






0.521058 290 


0.0134(8) 


-1/2 


0.016(6) 


0.02(4) 


-2 


10 




GE 


0.690 473 725 








0.22(2) 


-2 


16 






0.690 473 725 








0.18(4) 


-2.0(1) 


8 


Re 




0.690 5(6) 


0.13(3) 


-0.53(6) 








120 






0.690 473 725 


0.13(3) 


-0.53(3) 




0(7) 


-2 


120 




CE 


0.690 473 725 


0.10(3) 


-0.50(3) 


0.13(4) 






40 






0.690 473 725 


0.12(3) 


-0.52(3) 


0.05(9) 


0.9(5) 


-2 


20 






0.690 473 725 


0.102(2) 


-1/2 


0.10(2) 


0.5(3) 


-2 


24 
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FIG. 2: (Color online). Plot of /S.Rx = Rx — Oq — iJiL^'^''" versus L^^'^''" for the canonical wrapping 
probability R^. The straight line represents the leading finite-size dependence of IS.Rx, i.e. B2L?'^'"^'^. 
Parameter Oq and yean are fixed at the exact value and —1/2 respectively, and the amplitudes B\, 
B2 are obtained from the fit. 

probabilities Rb, Rx, and Re); further, for Rx the amplitude of Bi is very small. As 
expected, including terms with Bi and/or B2 only yields messy information and does not 
improve the quality of the fits. In the canonical ensemble, we also fitted the data by Eq. ( !22|l 
with a single correction term BiL^"'^'^ {B2 = Bi = 0). As shown in Tab. [H the data up 
to sufficiently large L^i^ (except Rx) have to be discarded for a reasonable residual x^- 
Nevertheless, we find that (i), the correction term has an exponent close to —1/2, and (ii), 
the estimates of Oq agree well with the exact values. A better description of the R data 
can be obtained by including correction terms with B^ and/or B2; for simplicity the exact 
values of Oq are used. As a result, the estimates of yean become more accurate, and are in 
good agreement with the prediction y^an = —1/2. We note that the correction coefficient 
Bi takes different values in the grand-canonical and the canonical ensemble. This is also in 
agreement with the theoretical expectation. As predicted in Sec. Ill, in addition to those 
with exponent n(2yt — d) {n = 1,2, . . .), correction terms with exponent —d — ip can also 
exist in the canonical ensemble. For the wrapping probabilities, one has = and thus 
such a correction term has the same exponent as BiL^'. Finally, in order to illustrate the 
term with L^^'^^"(= L~^), we plot ARx = Rx — Oq — -BiL^^'^" versus L^^'^''" as in Fig. [2l 
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TABLE II: Fit results for Si, S2, and the universal ratios Qs and Q 



Fit using O = L'^(Oo + BiL^-" + Br-Ly- + A^^O 





Oo 




l/can 


Bj- 






Si 


GE 


0.98448(3) 


- 


- 


1.7(13) 


-2(2) 


32 




CE 


0.9843(4) 


0.066(14) 


-0.50(5) 


6(5) 


-7(6) 


24 


S2 


GE 


1.035 97(4) 


- 


- 


1.7(7) 


-3.0(14) 


32 




CE 


1.035 8(5) 


0.05(3) 


-0.51(12) 


5(3) 


-7(5) 


24 


Qs 


GE 


0.96017(1) 






1.4(8) 


-1.8(12) 


40 




CE 


0.960 0(3) 


0.072 (8) 


-0.48(4) 






200 






0.960 03(13) 


0.075(5) 


-0.49(2) 


-0.3(16) 


0(2) 


24 


Qm 


GE 


0.870 56(2) 






-1.1(3) 


1.92(6) 


32 




CE 


0.870 4(2) 


0.138(7) 


-0.491 (13) 


-1.1(6) 


1.3(9) 


16 



In order to demonstrate that the correction exponent for the grand-canonical wrapping 
probabihties is indeed i/i = —2, we also fitted the data by Eq. fl22|) with i/i being free, Oq 
being fixed at the exact values and Bi = B2 = 0. The fit results are also included in Table [H 
The result is -2.02(4), -2.06(8), -2.4(6) and -2.0(1), for Ri, Rb, R^ and R^, respectively. 



These results are in compatible with the existing numerical result t/i = —2.003 (5) |18| for 
Ri, and agree well with the exact value yi = —2. It is interesting to mention that although 
correction of form L~'^{l + AlogL) [A is a constant) has been observed for other observables 



in the square-lattice bond percolation [2J], there is no evidence for such a logarithmic factor 
in the finite-size corrections for the wrapping probabilities. 

For magnetic quantities, we consider the largest-cluster size 5*1 and the second cluster-size 
moment S2, as well as the dimensionless ratios. We fitted the data by 

o = L'^iOo + Biiy-^-^ + Briy^- + Biiy'), (23) 

with Ui = —2 and ip being fixed at the exactly known value for the two-dimensional 
percolation. The term with L^'' is the correction from the regular background, with 
i/r = —yh = —91/48 for 5*1 and Qs, and yr = d — 2yh = —43/24 for S2 and Qm- The 
results are shown in Table HTl Again, we find that (i), the values of Oq remain unchanged 
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TABLE III: Fit results for the cluster-number density pk as determined for the bond (BP) and site 



(SP) percolation problem. Some fits make use of the exact results p\°q'^ = 0.098 076 211 
hg = 0.883 576 308 Q- 



2M and 



Fit using pk = pkfl + L~^{b + Biiy--- + Biiy^ 





Pk,0 


b 




yean 


Bi 






GE 


0.098 076 22(7) 


0.883 6(6) 


- 


- 


0.18(3) 


6 


BP 




0.098 076 211 


0.883 576 308 






0.178(14) 


6 






0.098 076 21(4) 


0.23(2) 


0.57(5) 


-0.47(6) 


0.0(4) 


8 




CE 


0.098 076 211 


0.226 (8) 


0.564(12) 


-0.47(3) 




10 






0.098 076 211 


0.233 (4) 


0.60(2) 


-1/2 




24 




GE 


0.027598 03(2) 


0.883 5(2) 






-0.17(4) 


16 


SP 




0.027598 02(2) 


0.883 576 308 






-0.22(3) 


20 




CE 


0.027598 00(4) 


0.41(2) 


0.64(7) 


-0.46(7) 




24 






0.02759799(4) 


0.417(2) 


0.687(11) 


-1/2 




24 



in the canonical ensemble, irrespective of whether they are universal (for Qs and Qm) or 
non-universal (for 5*1 and S2), and (ii), new correction terms with an exponent i/can = — 1/2 
are introduced. The estimate Qm,o = 0.870 56(2) in the grand-canonical ensemble agrees 
with the existing result 0.870 53(2) Q- 

B. Change of the universal excess cluster number 

For energy-like quantities, we consider the cluster- number density pk, and analyze the 
data by 

Pk = Pk,o + L-\h + B^iy- + B^Ly^), (24) 

De exactly obtained as p\^q'^ = (3-\/3 — 



2l|. In the grand-canonical ensemble. 



with Ui = —2 being fixed. The background term can 
5)/2 for the critical square-lattice bond percolation 
the excess cluster number b is known to be universal and the value has been obtained as 
bg = 0.883 576 308 19|, with subscript g for the grand-canonical ensemble. 



The fit results are given in Table lllli As in the above subsection, the correction term 
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CL 



0.09816 



0.09814 



0.09812 



0.0981 



0.09808 



grand-canonical 
canonical 




2e-05 



4e-05 
^-2 



6e-05 



8e-05 



FIG. 3: (Color online). Plot of the cluster-number density pk versus L ^ for the bond percolation 
model. The straight lines represent the leading finite-size dependence of pk as obtained from the 
fit. 

with Ucan = —1/2 is observed in the canonical ensemble, and the background contribution 
Pkfl remains unchanged. However, the excess cluster number is now be = 0.233(4), clearly 



different from the universal value bg = 0.883 576 308 [19|] in the grand-canonical ensemble. 
This indeed confirms the theoretical prediction for the case of ^ = —2 in Sec. Illli In other 
words, one expects that bg = bc + Bp'j^^, where B = p{l — j))/4 is a constant and p'^^ is the 
second derivative p'^ ^ of the regular part of pk with respect to the bond density. Since both 
B and p'l ^ are non-universal, the canonical-ensemble value b^ is no longer universal. As an 
illustration, we plot versus in Fig. [3l where the difference of bg and be is reflected by 
the different slopes of the data lines. 

In order to examine the "non-universal" nature of the excess cluster number be in the 
canonical ensemble, we also performed simulations of square-lattice site-percolation problem 



at the percolation threshold pf'^'^ = 0.592 746 02 ay ■ l27l|. In the grand-canonical ensemble, 
the simulations used 15 system sizes in range 4 < L < 512, and the results of the fits by 
Eq. ( l24l) are given in Tab. Illli The estimate of bg = 0.883 5(2) agrees well with the universal 



value 0.883 576 308 



19| 



In the canonical ensemble, the total number of occupied sites Ns{L) = p^^^'^L? is fixed. 
However, Ns{L) is not an integer, and thus the actual simulations were carried out for the 
total occupied-site number [Ns{L)] and [A^s(-^)] + 1, with [ ] for the floor integer. The 
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Monte Carlo results at Ns{L) were then obtained by linear interpolation. The simulations 
used 16 system sizes L in range 4 < L < 1024. The results of fits using Eq. fl2^ are 
also given in Tab. IIIII As expected, one has a correction term with exponent —1/2 and 
^(sitc)|_ Q 4]^7^2)] ^ bg. The "non-universal" property of is demonstrated by the fact 

^(sitc) _^ ^(bond)_ 

In addition, from our fits, one gets pf^^ = 0.027 598 00(5), which is in good agreement 



with the existing result 0.027598 1(3) 20|], and reduces the error margin significantly. 



We also study the FSS of the specific-heat-like quantity C2 for the bond percolation 
model and find that due to the limited precision, the data for L > 8 are well described by 
C2{L) = C2fi+cL~^ . As already mentioned in Section[IIl the excess fiuctuation c is universal. 
In the grand-canonical ensemble, the fit yields 6*2,0 = 0.039 446 (4) and Cg = 0.105 5(7); the 
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the latter agrees well 



former is consistent with the existing result 6*2,0 = 0.039 44(4) 
with the exactly known value c = 0.105 436 634 19|]. In the canonical ensemble, the result 
is C*2,o = 0.039 446(6) and Cc = —0.33(5). It is interesting to observe that not only the 
magnitude of the excess fiuctuation c is modified, but that also its sign has changed. 



C. Universality of scaling functions 

In Sec. Illlt we state that the scaling functions in Eqs. ( 1T9|) and ( l20l) are identical up to some 
correction terms, which refiects the universality of the functions. In order to demonstrate 
this, we carried out simulations near the critical point. We plot the grand-canonical 
versus L^* (p — p^) and the canonical versus L^' (p — Pc) as in Fig. HJ As expected, data 
points in both the ensembles nicely collapse to the same curve. 



V. DISCUSSION 



We derive the critical finite-size scaling behavior of percolation under the constraint that 
the total number of occupied bonds/sites is fixed, and confirm theoretical predictions for 
the two-dimensional percolation by means of Monte Carlo simulation in two dimensions. In 
particular, it is found that with the constraint, new finite-size corrections with exponent 
n{2yt — d) (n = 1, 2, . . .) are induced and the excess cluster number becomes non-universal. 
We note that our theory in Sec. IIIII can be used to explain the observed correction expo- 
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FIG. 4: (Color online). Plot of R X versus t for bond percolation. Parameter t represents L^'^{p—pc) 
in the grand-canonical ensemble (GE), and L^* (p— Pc) canonical ensemble (CE). The excellent 

collapse demonstrates the universality of scaling functions near the critical point. 

nent ^ —0.5 for the canonical wrapping probabilities in simulating the two-dimensional 



percolation by the Newman-Ziff algorithm 17|, ll8|. The predictions should be valid in any 



dimension with d > 2. We believe that this work provides an additional useful reference 
for percolation — a pedagogical system in the field of statistical mechanics. Furthermore, 
our work can help to understand the critical finite-size-scaling properties of other statistical 
systems in the canonical ensemble with a fixed total number of particles, which is the usual 
situation during experiments. 
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